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Beschreibung 

[0001] Die Erfindung betrifft ein Verfahren zur Bestimmung voo Frequenzen und Amplituden eines zumindest zwei 
Frequenzen enthaltenden Mehrtonsignals. 
5 [0002] Weiters belrifFt die Erfindung eine Vorrichtung zur Bestimmung von Frequenzen und Amplituden eines zumin- 
dest zwei Frequenzen enthaltenden Mehrtonsignals. - ■ ■ 
[0003] Mehrtonsignale, insbesondere Zweitonsignale, werden haufig eingesetzt, um Signalisierungsinformationen zu 
Uberlragen, beispielsweise bei der (Jbertragung von Wahlinformationen bei Telefonen im Sprachband, bei Zwischenarnt- 
signalisierungen, bei Femsteuerungsanwendungen (Anrufbeantworter . . .), etc. Aufgabe der Erfindung ist es. einen Weg 
10 zu schaffen, Frequenzen und Amplituden eines Mehrtonsignals, insbesondere eines Zweitonsignals zu bestimmen, um 
anhand dieser Signalparameter die Signalisierungsinfonnation zu erhalten und einer weiteren Verarbeitung zuganglich 
zu machen. 

[0004] Diese Aufgabe wird mit einem Verfahren der eingangs genannten Art erfindungsgemafi dadurch gelost, dass an- 
hand des Mehrtonsignals eine der Anzahl von Frequenzen entsprechende Anzahl von Filterkoeffizienten von Notchfil- 

15 termodellen berechnet wird, wobei die Filterkoeffizienten so berechnet werden, dass bei Durchgang des Mehrtonsignals 
durch eine durch die Filterkoeffizienten definierte Kombination von Notchfiltermodellen die Ausgangsleistung minimal 
ist. Mit dem erfindungsgemaBen Verfahren lassen sich sowohl die Frequenzen als auch die Amplituden eines Zweitonsi- 
gnals sowie die Energie des dem Zweitonsignal iiberlagerten Storsignales auf einfache Weise eindeutig bestimmen. 
[0005] Vorteilhafterweise wird das Zweitonsignal, wenn es als analoges Signal vorliegt in eine digitalisierte Form in 

20 Form eines Datenvektors (Xn) gebracht, aus welchem eine Matrix (D) der Form 



25 



Xj^ 4" Xq JC| 



gebildet wird. 

[0006] In einem nachsten Schritt konnen aus den noch tinbekannten Filterkoeffizienten der einzelnen Notchfiltenno- 
delle zumindest zwei Filterkoeffizienten der Kombination der Note hfiltermodelle gebildet werden, und aus den Filterko- 
30 effizienten der Kombination von Notchfiltermodellen kann ein Vektor der Form b = [1, bi, b2]^ gebildet werden. Die 
Komponenten bl und b2 dieses Vektors stellen die noch unbekannten Filterkoeffizienten der Kombination von Notchfil- 
termodellen dar. 

[0007] Liegt das Mehrtonsignal als Zweitonsignal vor so konnen fiir zumindest zwei Filterkoeffizienten der Kombina- 
tion der Notchfiltermodelle die Beziehungen bi = ai + 32 und b2 = aia2 + 2 gelten. 

35 [0008] Zur Bestimmung der Filterkoeffizienten kann eine quadratische Form der Gestalt b^D^Db minimiert werden. 
[0009] Die Frequenzen des Mehrtonsignals konnen nach den Formeln ©i = arccos(-ai/2), (O2 = arccos(-a2/2) und = 
arccos(-aii/2) berechnet werden, wobei ai, a2 und an die Filterkoeffizienten der einzelnen Notchfiltermodelle sind. 
[0010] Ist das Mehrtonsignal ein Zweitonsignal so hangen die Filterkoeffizienten der beiden einfachen Notchfiltermo- 
delle mit den Filterkoeffizienten der Kombination der Notchfiltermodelle bzw. des Doppelnotchfiltennodells nach fol- 

40 genden Formeln zusammen: 



45 -bjl^^b^ IA^2-b^ . 

[0011] Zur Bestimmung der Amplituden kann der Datenvektor (xj in folgende Form gebracht werden: 
Xn = ci cos(ncOi) + Si sin(na)i) + C2 cos(na>2) + S2 sin(nC02) + Eq, wobei n = e N gilt, 

und ci, si, C2, S2 Amplituden werte und die En kleine Fehlerwerte sind, wobei die Amplitudenwerte Ci, si, C2, S2 so be- 
50 stimmt werden, dass die Summe der Quadrate der En minimal wird und der Wert dieser Fehlerquadratsumme die Energie 
der iiberlagerten St5rung angibt. Uber die Frequenzen gehen die Filterkoeffizienten in die Form Xn = Ci cos(no)i) + s\ 
sin(nO>i) + C2 cos(nO>2) + S2 sin(na>2) + en, des Datenvektors ein. D. h. bei Kenntnis der Filterkoeffizienten konnen sowohl 
die Frequenzen als auch die Amplituden des Mehrtonsignals sowie die Energie der dem Mehrtonsignal iiberlagerten S to- 
rung berechnet werden. 

55 [0012] Vorteilhafterweise werden die Amplituden des Zweitonsignals nach folgenden Formeln berechnet: 

a, =^c/ -vs^ und 



60 a^zz^Ja^s^^ . 

[0013] Zur Durchfuhrung des erfindungsgemaBen Verfahrens eignel sich insbesondere eine Vorrichtung der eingangs 
genannten Art, welche dazu eingerichtet ist, anhand des Mehrtonsignals eine der Anzahl von Frequenzen entsprechende 
Anzahl von Filterkoeffizienten von Notchfiltermodellen zu berechnen, wobei sie dazu eingerichtet ist, die Filterkoeffi- 
65 zienten so zu berechnen, dass bei Durchgang des Mehrtonsignals durch eine durch die Filterkoeffizienten definierte 
Kombination von Notchfiltermodellen die Ausgangsleistung minimal ist. 

[0014] In einer bevorzugten Ausflihrungsform der Erfindung ist die Vorrichtung dazu eingerichtet, ein als analoges Si- 
gnal vorliegendes Zweitonsignal in eine digitalisierte Form in Form eines Datenvektors zu bringen, und aus dem Daten- 
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vektor eine Matrix der Form 



zu bilden. 

[0015] In einer bevorzugten Ausfuhrungsform der Erfindung kann die Vorrichtung dazu eingerichtet sein, aus den noch 
unbekannten Filterkoeffizienten der einzelnen Notchfiltermodelle zumindest zwei FilterkoefRzienten der Kombinalion lo 
der Notchfiltermodelle zu bilden, und aus den Filterkoeffizienten der Kombinalion einen Vektor der Form b = [1, bi, b2]^ 
zu bilden. wobei bl und b2 die noch unbekannten Filterkoeffizienten der Notchfilterkombination bedeuten. 
[0016] Liegt das Mehrtonsignal als Zweitonsignal vor so gelten fur die zwei Filterkoeffizienten der Kombinalion der 
Notchfiltermodelle die Beziehungen b^ = + ^2 ^"d b2 = aia2 + 2. 

[0017] Weitere Vorteile lassen sich dadurch erzielen, dass die Vorrichtung dazu eingerichtet ist, eine quadratische Form 15 
der Gestalt b^D^Db zur Bestimmung der Filterkoeffizienten der Kombinalion der Notchfiltermodelle zu minimieren. 
[0018] Welters kann die Vorrichtung dazu eingerichtet sein, die Frequenzen des Mehrtonsignals nach den Formeln ©1 
= arccos(-ai/2), a>2 = arccos(-a2/2) und 0)„ = arccos(-an/2) zu berechnen, wobei ai, ai und an die Koeffizienten der ein- 
zelnen Notchfiltermodelle sind. 

[0019] Wenn das Mehrtonsignal ein Zweitonsignal ist, kann die Vorrichtung dazu eingerichtet sein, die Filterkoeffi- 20 
zienten der einzelnen Notchfiltermodelle nach folgenden Formeln aus den Filterkoeffizienten der Kombinalion der 
Notchfiltermodelle zu berechnen: 



a,=6,/2-V6,V4 + 2-A2 , 



^2 =V^ 



25 



[0020] Eine bevorzugte Ausfuhrungsform der Erfindung sieht vor, dass die Vorrichtung dazu eingerichtet ist, die Stor- 
signal-Energie der dem Mehrtonsignal iiberlagerten Storung, aus den Filterkoeffizienten zu berechnen. 30 
[0021] Zur Bestimmung der Amplituden kann die Vorrichtung dazu eingerichtet sein, den Datenvektor in die Form Xn 
= ci cos(nC0i) + Si sin(nO>i) + C2 cos(n0)2) + S2 sin(nC02) + en zu bringen, wobei n = e N gilt, und ci, Si, C2, S2 Amplitu- 
denwerte und die Cn kleine Fehlerwerte sind, wobei die Vorrichtung welters dazu eingerichtet sein kann, die Amplituden- 
werte Ci, Si, C2, S2 so zu bestinunen, dass die Summe der Quadrate der minimal ist und der Wert dieser Fehlerquadrat- 
summe die Energie der dem Mehrtonsignal iiberlagerten Storung angibt. 35 
[0022] Daruber hinaus kann die Vorrichtung dazu eingerichtet sein, die Amplituden nach folgenden Formeln aus den 
Amplitudenwerten ci, si, C2, S2 zu berechnen: 



40 



[0023] Die Erfindung samt weiterer Vorteile wird im folgenden anhand einiger nicht einschrankender Ausfuhrungsbei- 
spiele, die in der Zeichnung veranschaulicht sind, schematisch daigestellt. In dieser zeigen: • 45 

[0024] Fig. 1 zwei Notchfiltermodelle (mit den Koeffizienten al und a2, welche die Sperrfrequenzen bestimmen) und 
[0025] Fig, 2 eine erfindungsgemaBe Vorrichtung zur Berechnung von Filterkoeffizienten der Notchfiltermodelle aus 
Fig. 1. 

[0026] GemaB Fig. 1 geht man zur Bestimmung der Amplituden ai, 32 und Frequenzen coi, CO2 eines Mehrtonsignals 
MSI von der Annahme aus, das Mehrtonsignal wiirde durch eine Hintereinanderschaltung von zwei Notchfiltermodeilen 50 
(Doppelnotchfilter) geschickt, wobei man aus dem gegebenen Eingangssignal die Koeffizienten al und a2 dieser beiden 
Notchfiltermodelle so ermittelt, dass bei der (nur als Hilfs annahme dienenden) Filterung des Zweitonsignales mit den 
hintereinandergeschalteten Notchfiltem die Ausgangsleistung A^ im eingeschwungenen Zustand minimal wird. Anhand 
der berechneten Filterkoeffizienten al und a2 kbnnen dann die Signalfrequenzen, die Amplituden sowie die Eneigie der 
dem Zweitonsignal iiberlagerten Storung ermittelt werden. 55 
[0027] GemaB Fig. 2 kann eine erfindungsgemaBe Vorrichtung VOR, falls das Mehrtonsignal in anaLoger Form vor- 
liegt, zunachsteine Tiefpassfilterung mittels eines Antialiasing-Filters ATFdurchfuhren und hiemach kann das Mehrton- 
signal einem A nalog-Digitai-Kon verier ADC zugefuhrt werden. 

[0028] Der Analog-Digital- Konverter ADC liefert Abtastwerte des Mehrtonsignals mit einer Abtastfrequenz, welche 
bekannterweise mindeslens doppelt so hoch sein muss wie die hochste zu verarbeitende Frequenzkomponente des Mehr- 60 
tonsignals. Beispielsweise betragt bei DTMF-Zcichen die hochste zu verarbeitende Signalfrequenzkomponente 
1633 Hz, sodass hier die Abtastfrequenz zu 4 kHz gewahlt werden kann. Der Analog-Digital- Wandler ADC iibermittelt 
seine Abtastwerte an eine Auswerteeinheit ASW zur Bestimmung der Amplituden aj, §2 und Frequenzen CDi, 0>2 des 
Mehrtonsignals MSI sowie zur Bestimmung der Enetgie N der dem Zweitonsignal iiberlagerten Storung. 
[0029] Die erfindungsgemaBe Auswerteeinheit ASW kann, beispielsweise ein digitaler Signalprozessor bzw. ein Mi- 65 
kroprozessor sein, der in einer zugehorigen Assemblersprache bzw. Hochsprache, wie Madab, Mathematica, etc, pro- 
grammiert sein kann. Die Abtastwerte konnen in einem Pufferspeicher PUS zwischengespeichert werden. Hierzu kann 
der Pufferspeicher PUS fUr 2 ♦ N Werte zykiisch beschrieben werden, wobei N die Blocklange eines Datenblockes ist, der 



3 
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auf das Vorhandensein eines Zweitonsignals ZSI untersucht wird! 

[0030] Solange die Auswerteeinhcit ASW eine Halfte des Pufferspeichers mit aktuell anfallenden Abtastwerten be- 
schreibt, konnen die Daten der andercn Pufferhalfte, wie im folgenden beschrieben, zur Bildung der in den folgenden 
Verarbeitungsschritten notigen GroBen verwendet warden (Ping-Pong-PufFer). Dadurch fallt mit der Taktfrequenz F/N 
5 jeweils ein neuer Datenvektor Xn an, der nun analysiert werden kann, wobei F die Abtastfrequenz ist und N beispiels- 

wdse zwischen__ca. 10 und 100-Werten liegt. - - - - - - — — ~ - 

[0031] Das erfindungsgemaBe Verfahren basiert auf der erweiterten Methode von Prony zur Schatzung der Parameter 
von im allgemeinen komplexen Exponendalfunktionen aus voigegebenen Datenwerten. 
[0032] Ein einfaches Notchfilter mit dem Filterkoeffizienten ai der Form 

10 

Hi(z) = 1 + aiz-^ + z-^ (1) 

kann eine Sinusschwingung sperren (j bedeutet im folgenden die positive Wurzel aus -1): 
15 1 + aie-J'*'^ + e-^j«^ = e-J"i (e*"! + aj + e-j"^) = e-j«i (ai + 2cosa>i) = 0 (2) 
wobei aus Gleichung (2) folgt 
ai = -2coso>i (3) 

20 

[0033] Die Kombination H(z) zweier Notchfilter Hi(z) und H2(z), 
H(z) = Hi(z)H2(z) (4) 
25 wobei 

H2(z)=l+a2Z-Uz-^ (5) 

gilt, und a2 der Filterkoeffizient des zweiten Notchfilters ist, kann ein Zweitonsignal ZSI sperren. 
30 [0034] Eine Kombination von Notchfiltem mit einer Anzahl von n Filtem, wobei gilt n > 0, n € IN, kann eine Anzahl 
von n Frequenzen sperren. 

[0035] Im folgenden soil das erfindungsgemalBe Verfahren anhand eines Zweitonsignals erlautert werden. Prinzipiell 
ist das erfindungsgemaBe Verfahren fiir Mehrlonsignale mit beliebig vielen Frequenzen geeignet, wobei jedoch der Re- 
chenaufwand mit zunehmender Anzahl der in dem Mehrtonsignal enthaltenen Frequenzen steigt. 
35 [0036] Mit (5) wird (4) zu 

H(z) = 1 + (ai + a2)z'^ + (a; 32 + 2)z"^ + (ai + a2)z~"^ + z"^ (6), 

wobei die zweite Sperrfrequenz durch 

40 

a2 = -2cos(02 (7) 

gegeben ist Mit den Hilfskoeffizienten 
45 bi = ai + a2 (8) 
■ b2 = aia2 + 2 (9) 

wird die Obertragungsfiinktion des Doppelnotchfilters (6) zu 

50 

H(z) = 1 + biz-^^ + b2Z-^ + biz'^ + (10). 

Zur Bestimmung der beiden Frequenzen des Zweitonsignals ZSI werden die beiden Hilfskoeflizienten bi und b2 in (10) 
so gew^lt, dass das durch seine Abtastwerte gegebene Zweitonsignal 2^1 

55 

x„,n = O...N-l (11) 

beim Durchgang durch die Kombination H(z) der Notchfilter Hi(z), H2(z) maximal gedampft wird. 
[0037] Die Auswerteeinheit ASW ist dazu eingerichtet, aus den Datenwerten Xn zunachst die Hilfskoeffizienten b 1 , b2 
60 des Doppelnotchfilters, dann die Filterkoeffizienten al, a2 der einzelnen Notchfilter, dann die Frequenzen ©i, 0)2, dann 
die Amplituden ai, %i sowie die Energie N der dem Zweitonsignal uberlagerten Storung zu berechnen, 
[0038] Zunachst wird das Signal am Ausgang des Modells des Doppelnotchfilters, das nur als Hilfskonstruktion zur 
Besdmmung der Filterkoeffizienten dient, in Abhangigkeit der Hilfskoeffizienten bl, b2 dargestellt. 

65 
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[0039] Fiir den Fall N = 6, d. h. die minimal notwendige Anzahl von Datenwerten zur Bestimmung der sechs Parame- 20 
ter des Zweitonsignals (2 Frequenzen, 2 Amplituden, 2 Nullphasenwinkel) konnen aus (12) bi und b2 bestimmt werden, 
in dem die beiden linearen Gleichungen fiir 



A4 = 0 (13) 
A5 = 0 (14) 



25 



geldst werden. Zur Verminderung des Einflusses eventuell vorhandener Rauschsignale in den Datenwerten ist es jedoch 
vorteilhaft N grofier als 6 zu wahlen. 

[0040] Die Koeffizienten bi und b2 werden so bestimmt, dass die Filterausgangsleistung fur den eingeschwungenen 30 
Zustand minimal wird, d. h. 



^A^^^vmn (15) 



[0041] Mit der Matrix 



(16) 



und dem Vektor 

b = [l,bi,b2f (17) 

fiihrt dies auf die Minimierung der quadratischen Form 

b'^D'^Db-*^min (18) 

mit der symmetrischen Matrix 





^00 


«01 


«Q2" 






E = 






«12 




(19) 




/20 


^21 


*22. 







[0042] Da die Matrix E aufgrund ihrer Konstruktion symmetrisch ist, miissen nur die auf und z. B. iiber der Hauptdia- 
gonale liegenden Elemente berechnet werden. 

[0043] Die quadratische Form (18) nimmt ihr Minimum fiir die Hilfskoeffizienten 

bi = (eo2ei2 - eoie22)/(eiie22 - Qxz^) (20) 

und 

= (eoiei2 - eiieo2)/(eiie22 - e^^) (21) 
an. Wenn der Ausdruck im Nenner (ene22 - ^ii^) dem Betrag nach sehr klein wird oder verschwindet, ist dies ein Zeichen 
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dafur, dass in dem analysierten Datenvektor Xn kein ausreichend storungsfreies Zweitonsignal ZSI enthalten ist Dies ist 
z. B. bei einem einzigen Sinussignal der Fall (keine zweite Sinuskomponente vorhanden). Ist hingegen der Nenner dem 
Betrag nach hinreichend groB, so werden die Hilfskoeffizienten bi und b2 berechnet. 

[0044] Bestehl das betxachtete Datensignal aus zwei reinen Sinussen (d. h. es ist kein Storsignal uberlagert), wird das 
5 Minimum zu Null, da dann das Filter H(z) beide Frequenzen exakt spent. Lost man die Gleichungen (8) und (9) nach den 
Filterkoeffizienten ai und ai der einzelnen Notchfilter auf, so erhalt man: - 

a, =fr,/2-VA,'/4+2-Aj (22) 

10 a, =ft,/2 + ^A,'/4 + 2-A2 (23), 

[0045] In der Auswerteeinheit ASW kann gepruft werden, ob der Ausdruck unter der Wurzel in (22) und (23) positiv 
ist. Ist dies nicht der Fall, so kann kein Zweitonsignal ZSI von der Auswerteeinheit detektiert werden. Weiters kann die 
Auswerteeinheit ASW dazu eingerichlet sein, zu iiberprufen, ob die beiden Werte ai und 32 im Intervall [-2, 2] liegen, 
15 was zur Umkehrung der Gleichungen (3) und (7) erforderlich ist. Ist diese Bedingung nicht erfullt, so wird von der Aus- 
werteeinheit ASW kein Zweitonsignal ZSI detektiert. Mit (3) und (7) konnen die Werte fiir die beiden Signaifrequenzen 



= arccos(-ai/2) (24) 

20 

und 

(i>2 = arccos(-a2/2) (25) 



25 berechnet werden. 

[0046] Allgemein lasst sich sagen, dass sich im Fall einer Anzahl von n Frequenzen in dem Mehrtonsignal die n-te Fre- 
quenz (On aus dem n-ten Filterkoeffizienten an des die n-te Frequenz <»a sperrenden einfachen Notchfiltermodelles nach 



30 



35 



<Dn = arccos(-an/2) (25a) 
berechnen lasst. 

[0047] Die zugehorigen Frequenzen analoger Sinuskomponenten ergeben sich zu 
/i = F • coi/(2 • 71), /2 = F . ©2/(2 • 71) und /n = F • ay(2 • n). 



[0048] Da die Berechnung der inversen Cosinusfunktion numerisch aufwendig ist, wird sic nur dann durchgefuhrt, 
wenn es unbedingt notwendig ist. FQr die Bestimmung der Amplituden und der StOrsignal-Energie kann man darauf ver- 
zichten. Wenn nicht die numerischen Werte der ermittelten Frequenzen als Endergebnis benotigt werden, sondern nur das 
Ergebnis eines Vergleichs der ermittelten Frequenzen mit ihren Sollwerten, so braucht man die inversen Cosinusfiinktio- 
40 nen nicht zu berechnen, da die Vergleichsoperationen auch zwischen den Kosinuswerten der gemessenen Frequenzen 
und den Kosinuswerten der entsprechenden Schwellen durchgefuhrt werden konnen. 

[0049] Zur Bestimmung von Amplitudenwerten Ci, Si, C2, S2 kann fiir die gegebenen Datenwerte folgender Ansatz 
gemacht werden: 

45 Xn = ci cos(n©i) + Si sin(n(Di) + C2 cos(nC02) + Sz sin (niOi) + 8n, n = 0, . . ., N - 1 (26), 

wobei die Cn kleine Fehlerwerte sind. Die Amplitudenwerte ci, Sj, C2, S2 werden in der Auswerteeinheit ASW so be- 
stimmt, dass die Summe der Quadrate der minimal wird. 
[0050] IVIit 



50 



55 



X = [xo, xi, . . ,xn-i]'^ (27), 

cos(Oa>,) sin(OG;|) 
cos(lca,) sin(l£0,) 



cos(0^ 2) 
cos(lct?2) 



sin(0^2) 



cos((N-l)a),) sin(iN-\)co,) cos((JV -!)£»,) sin((JV-l)6>,) 

60 C = [Ci,Si,C2,S2]'^ (29) 

und 

e=[eo,ei....eN-i]'^ (30) 



,i = 0..JV-l,y = 0...3 (28; 



65 



kann (26) in Matrixschreibweise als 



x-Bc = e (31) 
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geschrieben werden. 

[0051] Die Elemente der obersten Zeile von B sind 1 fur j = 0, 2 (da cos(O) = 1) und 0 fur j = 1, 3 (da sin(O) = 0). Die 
Elemente der zweiten Zeile sind ebenfalls bekannt, da sie fur j = 0, 2 mit (3) und (7) sofort aus den schon bekannten Ko- 
effizienten_ajjjnda2jefunden werden kdnnen und sich die zugehorigen Sinuswerte.bei j = 1,3 sich durch die Beziehung 
sin(a») =yi - cos^(0))_berechnen lassen.-Die restlichen-Elemente der Spalten von B lassen sich mit Hilfe der Relcursions- 
gleicliung 

x(n) - 2cos((0)x(n - 1) + x(n - 2) = 0 
berechnen, welche fiir trigonometrische Folgen der Bauart 
x(n) = cos(n • (0) 

bzw. 15 
x(n) = sin(n • co) 

erfiillt ist, wobei die ersten beiden, nun schon bekannten Elemente jeder Spalte von B als Startwerte der Rekursion be- 

niitzt werden konnen. 20 
[0052] DieBedingung 

e'^e min (32) 

fuhrt auf das'lineare Gleichungssystem LGS: 25 
Gc = g (33) 

mil der auf Grund ihrer Konstruktion synunetrischen Matrix G und dem Vektor g 

30 

G = B'^B (33a) 
g = B'^x (33b) 



welches mit einem direkten oder einem iterativen Verfahren, beispielsweise dem Gauss'schen Verfahren, nach dem Ko- 35 
effizientenvektor c aufgelost werden kann. Aus den Komponenten von c erhalt man die endgiiUigen Amplitudenwerte 

3, =Vc,'+V (34) 

und 40 
«2 =V^^+?~ (35) . 

[0053] Die nachfolgende Bewertung der Amplituden Si, §2 kann aber auch mit den Quadraten der Amplituden ai, §2 
durchgefuhrt werden, weshalb in einer bevorzugten Ausfiihrungsform der Erfindung zugunsten eines verringerten Re- 
chenaufwandes in der Auswerteeinheit ASW auf die Berechnung der Wurzel in (34) und (35) verzichtet wird. 
[0054] Dariiber hinaus kann zur Bewertung der Amplituden Si, S2 iiberpriift werden, ob sie innerhalb geforderter Tole- 

ranzgrenzen liegen. 

[0055] Welters kann die Auswerteeinheit ASW dazu eingerichtet sein, die in dem analysierten Datenvektor Xn enthal- 
tene Storenergie N nach Gig. (36) zu berechnen, wobei die Storenergie N der Summe der Fehlerquadrate beim Least- 
Squares-Fii zur Berechnung der Signalamplituden entspricht. 

N = x'^x-g'^c (36) 

[0056] Anhand der Storenergie N kann das Signal-Rausch-Verhaltnis zur Signalbewertung beslirnmt werden. Sind alle 
Toleranzanforderungen fur den aktuell analysierten Datenvektor erfullt, so kann auch noch weilers uberpruft werden, 
ob die Toleranzanforderungen fiir eine vorgebbare Anzahl unmitlelbar vorangegangener Datenvektoren Xo erfUllt sind. 
Diese "Qberpriifung kann in der Auswerteeinheit ASW beispielsweise in Form einer sogenannten "State-Machine" imple- 
mentiert sein. Das Ergebnis ERG der Uberpriifung ist entweder "ein durch ein Mehrtonsignal, beispielsweise ein Zwei- 
tonsignal ZSI, reprasentiertes Zeichen wurde erkannt" oder "kein Zeichen erkannt". 

[0057] Die berechneten Amplituden und Frequenzen konnen von der Auswerteeinheit ASW auch einer Anwendung 
ANW bzw. einem Anwendungsprogramm zur weiteren Verwendung zur Verfugung gestellt werden. 

Patentanspriiche 

65 

1. Verfahren zur Bestimmung von Frequenzen (Oi, 0)2) und Amplituden (ai, a2) eines zumindest zwei Frequenzen 
(0)i, 0)2) enthaltenden Mehrtonsignals (ZSI), dadurch gekcnnzeichnet, dass anhand des Mehrtonsignals eine der 
Anzahl von Frequenzen (cOi, CO2) entsprechende Anzahl von Filterkoeffizienten (ai, a^) von Notchfiltermodellen 



12/23/04, EAST Version: 2.0.1.4 



DE 101 31 251 A 1 



(Hi(z), H2(z)) berechnet wird, wobei die Filterkoeffizienten (ai, 82) so berechnet werden, dass bei Durchgang des 
Mehrtonsignals durch eine durch die Filterkoeffizienten (ai, 32) definierte Kombination (H(z)) von Notchfiltermo- 
dellen (Hi(z), Eiiz)) die Ausgangsleistung (A^) minimal ist 

2. Verfahren nach Anspruch 1, dadurch gekennzeichnet, dass das Mehrtonsignal (ZSI) als analoges Signal vorliegt 
5 und in eine digitalisierte Form in Form eines Datenvektors (Xq) gebracht wird, aus welchem eine Matrix (D) der 

Form _ _ — - ------- - - 



gebildet wird. 

3. Verfahren nach Anspruch 1 oder 2, dadurch gekennzeichnet, dass aus den noch unbekannten Filterkoeffizienten 
15 (ai, BLi) der einzelnen Notchfiltermodelle (Hi(z), H2(z)) zumindest zwei Filterkoeffizienten (b^, b2) der Kojnbination 

(H(z)) der Notchfiltermodelle (Hi(z), H2(z)) gebildet werden, und aus den Filterkoeffizienten (bi, b2) ein Vektor (b) 
der Form b = [1, bi, b2]^ gebildet wird. 

4. Verfahren nach Anspruch 3, dadurch gekennzeichnet, dass fiir zwei Filterkoeffizienten (bi, b2) der Kombination 
(H(z)) der Notchfiltermodelle (Hi(z), H2(z)) die Beziehungen bi = ai + a2 und b2 = aia2 + 2 gelten. 

20 5. Verfahren nach den Anspriichen 3 und 4, dadurch gekennzeichnet, dass eine quadratische Form der Gestalt 

b^D^Db zur Bestimmung der Filterkoeffizienten (bl, b2) minimiert wird. 

6. Verfahren nach einem der Anspriiche 1 bis 5, dadurch gekennzeichnet, dass die Frequenzen (cOi, C02)des Mehr- 
tonsignals (ZSI) nach den Formeln a>i = arccos(-ai/2), CO2 = arccos(-a2/2) und co^ = arccos(-afl/2) berechnet werden, 
wobei ai, a2 und a,, die Koeffizienten der einzelnen Notchfiltennodelle sind. 
25 7. Verfahren nach Anspruch 6, dadurch gekennzeichnet, dass das Mehrtonsignal ein Zweitonsignal ist und Koeffi- 

zienten(ai, der einzelnen Notchfiltermodelle nach folgenden Formeln aus den Filterkoeffizienten (bi, b2) der 
Kombination (H(z)) der Notchfiltermodelle (Hi(z), H2(z)) berechnet werden: 

a, =6,/2-Ja'/4+2-62 , 

30 I I V I 

02 =6,/2 + V*i^/4 + 2-6, . 

8. Verfahren nach einem der Anspriiche 1 bis 7, dadurch gekennzeichnet, dass die Storsignal-Energie N der dem 
35 Mehrtonsignal iiberlagerten Storung, aus den Filterkoeffizienten (ai, 3.2) berechnet wird. 

9. Verfahren nach einem der Anspriiche 2 bis 8, dadurch gekennzeichnet, dass zur Bestimmung der Amplituden (ai, 
§2) der Datenvektor (Xn) in folgende Form gebracht wird: 

Xn = ci cos(n(Oi) + Si sin(n(Oi) + C2 cos(no>2) + S2 sin (n(02) + tn* wobei n = 0, . . N - 1 gilt, und ci, si, C2, S2 Am- 
plitudenwerte und die kleine Fehlerwerte sind. 
40 10. Verfahren nach Anspruch 9, dadurch gekennzeichnet, dass die Amplituden werte c^, si, C2, S2 so bestimmt wer- 

den, dass die Summe der Quadrate der etj minimal wird und der Wert dieser Fehlerquadratsumme die Energie der 
dem Zweitonsignal iiberlagerten Storung angibt. 

11 . Verfahren nach Anspruch 9 oder 10, dadurch gekennzeichnet, dass die Amplituden (§1, §2) nach folgenden For- 
meln berechnet werden: 

45 

12. Vorrichtung zur Bestimmung von Frequenzen (o)i, o>i) und Amplituden (ai, ^2) ^i^^s zumindest zwei Frequen- 
zen (coi, 02) enthaltenden Mehrtonsignals (ZSI), dadurch gekennzeichnet, dass sie dazu eingerichlet ist, anhand des 

50 Mehrtonsignals eine der Anzahl von Frequenzen ((Oi, 0)2) entsprechende Anzahl von Filterkoeffizienten (ai, a2) von 

Notchfiltermodellen (Hi(z), H2(z)) zu berechnen, wobei sie dazu eingerichlet ist, die Filterkoeffizienten (ai, a2) so 
zu berechnen, dass bei Durchgang des Mehrtonsignals durch eine durch die Filterkoeffizienten (ai, Hi) definierte 
Kombination (H(z)) von Notchfiltermodellen (Hi(z), HzCz)) die Ausgangsleistung (A^) minimal ist. 

13. Vorrichtung nach Anspruch 12, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, ein als analoges Signal 
55 vorliegendes Zweitonsignal (ZSI) in eine digitalisierte Form in Form eines Datenvektors (xn) zu bringen, und aus 

dem Datenvektor (Xn) eine Matrix (D) der Form 

X4 + X(j X3 + X, Xj 



.^W-l +^iV-5 ^A/-4 ; 

zu bilden. 

14. Vorrichtung nach Anspruch 12 oder 13, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, aus den unbe- 
65 kannten Filterkoeffizienten (aj, 32) der einzelnen Notchfiltermodelle (Hi(z), H2(z)) zumindest zwei Filterkoeffizien- 
ten (bl, hi) der Kombination (H(z)) der Notchfiltermodelle (Hi(z), H2(z)) zu bilden, und aus den Filterkoeffizienten 
(bl, b2) einen Vektor (b) der Fonn b = [1, bi, b2)^ zu bilden. 

15. Vorrichtung nach Anspruch 14, dadurch gekennzeichnet, dass fiir zwei Filterkoeffizienten (bi, b2) der Kombi- 



10 ^ = 



60 
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nation (H(z)) der Notchfiltermodelle (Hi(z), HaCz)) die Beziehungen bi = ai + a2 und b2 = aia2 + 2 gelten. 

16. Vorrichtung nach einem der Anspriiche 13 bis 15, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, eine 
quadratische Form der Gestalt b'^D^Db zur Bestimmung der Filterkoeffizienten (bi, b2) der Kombination (H(z)) der 
Notchfiltermodelle (Hi(z), H2(z)) zu minimieren. 

17. Vorrichtung nach einem der Anspriiche 12 bis 17, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, die 5 
Frequenzen ((Oi, (p2) des Mehrtonsignals (ZSD nach den Formeln-COi = arccos(-ai/2), CO2 =-arccos(=a2/2) und cOn = ' " 
arccos(-an/2) zu berechnen, wobei ai, a2 und an die Filterkoeffizienten der einzelnen Notchfiltermodelle (Hi(z), 
H2(z)) sind. 

18. Vorrichtung nach Anspruch 17. dadurch gekennzeichnet, dass das Mehrtonsignal ein Zweitonsignal ist und die 
Vorrichtung dazu eingerichtet ist, die Filterkoeffizienten (ai, ai) der einzelnen Notchfiltermodelle (Hi(z), H2(z)) 10 
nach folgenden Formeln aus den Filterkoeffizienten (bi, b2) der Kombination (H(z)) der Notchfiltennodelle (Hi(z), 
H2(z)) zu berechnen: 



19. Vorrichtung nach einem der Anspriiche 12 bis 18, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, die 
St£3rsignal-Energie (N) der dem Mehrtonsignal Oberlagerten Storung, aus den Filterkoeffizienten (aj, 82) zu berech- 
nen. 

20. Vorrichtung nach einem der Anspruche 13 bis 19, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, zur 
Bestimmung der Amplituden (§1, §2) den Datenvektor (Xq) in die Form Xn = Ci cos(na)i) + si sin(na)i) + C2 cos(nCD2) 20 
+ S2 sin(ncc>2) + En zu bringen, wobei n € N gilt, und ci, si, C2, S2 Amplitudenwerte und die 8n kleine Fehlerwerte 
sind. 

21. Vorrichtung nach Anspruch 18, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, die Amplitudenwerte 
(ci, Si, C2, S2) so zu bestimmen, dass die Summe der Quadrate der En minimal ist und anhand dieser Fehlerquadrat- 
summe die Storsignalenergie zu berechnen. 25 

22. Vorrichtung nach Anspruch 18 oder 19, dadurch gekennzeichnet, dass sie dazu eingerichtet ist, die Amplituden 
(ai, §2) nach folgenden Formeln aus den Amplitudenwerten (c^, Si, C2, S2) zu berechnen: 



15 
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30 
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